function res = BD_sin(omega,x,h)
res = (sin(omega*x)-sin(omega*(x-h)))/h;